Current noise in long diffusive SNS junctions in incoherent MAR regime. 
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Spectral density of current fluctuations at zero frequency is calculated for a long diffusive SNS 
junction with low-resistive interfaces. At low temperature, T ^ A, the subgap shot noise approaches 
linear voltage dependence, S = {2/3R){eV + 2A), which is the sum of the shot noise of the normal 
conductor and voltage independent excess noise. This result can also be interpreted as the 1/3- 
suppressed Poisson noise for the effective charge q = e(l-|-2A/eV^) transferred by incoherent multiple 
Andreev reflections (MAR). At higher temperatures, anomalies of the current noise develop at the 
gap subharmonics, eV — 2A/n. The crossover to the hot electron regime from the MAR regime is 
analyzed in the limit of small applied voltages. 
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During the last decade, considerable attention has been 
focused on the study of current fluctuations in meso- 
scopic systems, especially on shot noise which reflects 
the discrete nature of the electron charge and correla- 
tions between electrons. Whereas in normal ballistic sys- 
tems with tunnel barriers the spectral density S of the 
shot noise at zexo frequency approaches full Poissonian 
value Sp = 2elp mesoscopic diffusive wires shorter than 
the inelastic scattering length produce Poisson noise sup- 
pressed by a factor 1/3. The universality of this facUa 
was demonstrated within diffexent theoretical modelaa'El 
and confirmed experimentally.u The juelastic scattering 
could change the suppression factorJS'B which was used 
in experimentsQ as a probe of the electron relaxation. 

In hybrid normal-superconducting (NS) systems, the 
current transport at subgap voltages, eV < A, is associ- 
ated with the transfer of an elementary charge 2e, instead 
of e, due to Andreev reflection of quasiparticles from the 
NS interface, which converts two electrons in the normal 
metal to a Cooper pair in the superconductor. As shown 
theoreticallylZl and observed in experimentsjs such dou- 
bling of the elementary charge leads to doubling of the 
subgap shot noise in NS junctions. 

A more pronounced enhancement of the shot noise is 
expected in superconducting tunnel and SNS junctions 
whose subgap conductivity involves multiple Andreev re- 
flections (MAR). In this case, the effective transferred 
charge q = {Na + l)e increases along with the num- 
ber Na ^ 2A/eV of Andreev reflections of quasiparticle 
which overcomes the energy gap 2A by elementary steps 
eV. As a result, the shot noise spectral density S{V) at 
low voltages, eT^ ^ A, should greatly exceed the Poisso- 
nian value and approach a constant level, which can be 
estimated as S'(O) = 2qV/R = 4A/i? for ballistic junc- 
tions with normal resistance R, and as 

S{0) = AA/m (1) 

for diffusive junctions taking into account the 1/3- 



suppression factor. Alternatively, the enhanced noise 
may be interpreted as "thermal" noise generated by 
nonequilibrium quasiparticles within the whole subgap 
region \E\ < A, which is a characteristic property of the 
MAR regime. Q The experimental observation of multiply 
enhancecUshot noise was reported for NbN-based tunnel 
junctionst2l and several types of SNS junctionsJlJ 

Theoretical analysis of the shot noise in short junc- 
tions, with length d smaller than the coherence length 
fo = {hD/Ay^'^ in the superconductor, was done in Refs. 
n3 and n3. In such systems, the quantum coherence be- 
tween the electrons and retroreflected holes extends over 
the entire junction, which leads to the ac.Josephson ef- 
fect and to non-ohmic I-V characteristico Nevertheless, 
the estimate q ^ {Na + l)e for the effective transferred 
charge holds, although the zero-bias spectral density of 
the shot noise differs from Eq. Along with the dc 
current, the shot noise reveals subharmonic gap struc- 
ture (SGS), i.e., steps in S{V) or dS/dV at V = 2A/ne 
(n = 1, 2, . . .). A semiclassical model of the shot noise in 
ballistic point contacts was proposed in Rcf. |l^. 

In this paper we analyze the current noise in long 
diffusive SNS junctions, d ^ ^o- When the bias volt- 
age is much larger than the Thouless energy i^xh = 
fiD/d^ <^ A [D is the diffusion coefficient), then the 
size of the coherent proximity regions near the NS inter- 
faces, = {TlD / Ey/"^ , is much smaller than the junc- 
tion length d at all relevant energies E > eV . In this 
case, the Josephson effect is suppressed, and the sub- 
gap current transport can be-quantitatively described in 
terms of incoherent MAR,cHla similar to_the quasiclas- 
sical theory for ballistic systems (OTBKIlj). The shot 
noise in the incoherent MAR regime was analyzed in Ref. 
^ for a diffusive SNS junction with a tunnel barrier in- 
side the normal metal, assuming the barrier to dominate 
the junction resistance. Under these conditions, the shot 
noise is generated by tunneling electrons and cau. there- 
fore be calculated within the tunnel approach.LJ If the 



1 



resistance of the normal metal exceeds the resistance of 
possible tunnel barriers within the junction, shot noise 
emerges due to impurity scattering. In thiSj-oase, it can 
be calculated within a Langevin approachjla if we ne- 
glect the contribution of the small proximity regions in 
the vicinity of the interfaces. Following Ref. |l9|, in which 
the Langevin equation was applied to the current fluctu- 
ations in a diffusive NS junction, we derive an expression 
for the current noise spectral density in SNS junctions 
at zero frequency in terms of the nonequilibrium popu- 
lation numbers n'^'^{E,x) of electrons and holes within 
the normal metal, Q < x < d, 

S=^j'^fjE [n^ (1 - + (l - n^)] ■ (2) 
The electric current through the junction is given by 

(3) 



d f°° 
2eR J_^ 



The population numbers obey the diffusion equation. 
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with the inelastic collision term I^. At \E\ > A, the 
boundary populations nQ^d{E) = n{E,x) \x=o,d are local- 
equilibrium Fermi functions, Uq'^ {E)=np{E), n'^'^{E) — 
np{E ± eV) (we use the potential of the left electrode as 
the energy reference level). At subgap energies, \E\ < A, 
the boundary conditions should be modified in accor- 
dance with the mechanics of complete Andreev reflection 
which equalizes the electron and hole population num- 
bers at a given electrochemical potential and blOjCks the 
net probability current through the NS interface,t3 

n^E) = 4{E), nl{E - eV) = n'^^iE + eV), 
<(£;) + 4'{E) = 0, n'^JiE - eV) + n'i'{E + eV) = 0, (5) 

where tIq ^ are the boundary values of the electron and 
hole probability flows dn/dx. 

In the absence of inelastic scattering, the population 



numbers are linear functions of a;, n'^' (E, x 



\E)- 



xn'^'^'(E), which results in the recurrences for boundary 
populations and diffusive flows within the subgap region. 



nl'^{E - eV) - n^^iE + eV) = T^dn'^'^iE T eV), (6) 



e.h 



e.h/ 



{E ~ eV) ^ n/"{E + eV). 

According to Ref. |l5|, these recurrences are equivalent to 
the problem of "current" and "voltage" distribution in an 
equivalent network in energy space. In the present ap- 
proximation, which assumes the contribution of the prox- 
imity effect and the normal scattering at the interfaces 
to be negligibly small, this_aetwork consists of a series 
of resistances of unit valueo connected periodically, at 
the energies Ek = E + keV, with the distributed "volt- 
age source" npiE) (see Fig. ft). The "potentials" nk of 



the network nodes with even numbers k represent equal 
electron and hole populations nQ^{E + keV) at the left 
NS interface, whereas the potentials of the odd nodes de- 
scribe equal boundary populations n^/^i^E + keV =F eV) 
at the right interface. The "currents" Ik entering fc-th 
node are related to the probability currents n'{Ek) as 
h{E) = -dn'"{Ek-i) (odd k) and hiE) = dn'"{Ek) 
(even fc), and represent partial electric currents trans- 
ferred by the electrons and holes across the junction, 
obeying Ohm's law in energy space, Ik — nk-i — Uk- 
Within the gap, \Ek\ < A, i.e., at -iV_ < fc < iV+, 
N±{E) = Int[(A =F E)/eV] + 1 [Int(x) denoting integer 
part of x], the nodes are disconnected from the reser- 
voir due to complete Andreev reflection and therefore all 
currents flowing through the subgap nodes are equal. 

Due to periodicity of the network, the partial currents 
obey the relationship Ik{E) = Im[E+{k — m)eV], and the 
boundary population hq is related to the node potentials 
Uk as no{E + keV) — nk{E). This allows us to reduce 
the integration over energy in Eqs. (|^) and (^ to an 
elementary interval < E < eV, 



k— — oo 



eV 



k—~oo 



2nk{l - Uk) + -II 



(7) 
(8) 



1-1 'o li 



-N 




-A 







A 



FIG. 1. Equivalent MAR network in energy space in the 
limit of negligibly small normal reflection and proximity re- 
gions at the NS interfaces. 

The "potentials" of the nodes outside the gap, \Ek\ > 
A, are equal to local-equilibrium values of the Fermi func- 
tion, nk{E) = np{Ek) at fc > A^+, k < —N^. The partial 
currents flowing between these nodes. 



h = npiEk-i) - npiEk), k > N+, k < -N^ 



(9) 



are associated with thermally excited quasiparticles. The 
subgap currents may be calculated by Ohm's law for the 
series of + subgap resistors. 



Ik = 



n+ 



N+ + N_' 



-N- <k< N+, 



(10) 



where n±{E) — nFiE±N±)- From Eqs. (||) and (|T^) 
we obtain the current spectral density in Eq. (^ as 
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Sfc-oo -^fe = 1' which results in Ohm's law, V = IR, 
for the net electric current through the junction. This 
conclusion is closely related to our disregarding the prox- 
imity effect and the normal scattering at the interface. 
Actually, both of these factors lead to the appearance 
of SGS and excess or deficit currents in the I-V char- 
acteristic, with the magnitude increasing along with the 
interface barrier strength and the ratio / dS3 

The subgap populations can be found as the potentials 
of the nodes of the subgap "voltage divider" , 
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(11) 



By making use of Eqs. the net current noise 

can be expressed through the sum of the thermal noise 
of quasiparticles outside the gap. 



eV 



In 



npiA + eV) 



«f(A) 



coth 



eV 

2r 



(12) 



and the subgap noise. 



Sa = ^£ dE{N++N-) [/+_ + /_+ 
+2(/++ + /__)], faf3 = n4l-n0). 



(13) 



At low temperatures, T ^ A, the thermal noise Sy 
vanishes, and the total noise coincides with the subgap 
shot noise, which takes the form 



5* 



2 

3^ 



dE{N+ + N-) 



3i? 



(ey + 2A), (14) 



of 1/3-suppressed Poisson noise S = {2/3)qI for the ef- 
fective charge q = e{l + 2A/eV). At ^ ^ 0, the shot 
noise turns to a constant value 4A /3i? in Eq. , which is 
identical to the result of Ref. ^ and therefore seems to be 
universal for the incoherent MAR regime in long diffusive 
junctions and independent of the shot noise mechanism. 
At finite voltages, this quantity plays the role of the "ex- 
cess" noise, i.e. the voltage-independent addition to the 
shot noise of a normal metal at low temperatures [see 
Fig. §(a)]. Unlike short junctions, where the excess noise 
is proportional to the excess current in our system the 
excess current is small and has nothing to do with large 
excess noise. 

Results of numerical calculation of the noise at finite 
temperature are shown in Fig. ^. While the temperature 
increases, the noise approaches its value for normal metal 
structures,0 with additional Johnson-Nyquist noise com- 
ing from thermal excitations. In this case, the voltage- 
independent part of current noise may be qualitatively 
approximated by the Nyquist formula S{T) ~ AT* / R 
with the effective temperature T* ^ T + A(T)/3. The 
most remarkable phenomenon at nonzero temperature is 



the appearance of steps in the voltage dependence of the 
derivative dS/dV at the gap subharmonics eV = 2A/n 
[Fig. §(b)], which reflect discrete transitions between the 
quasiparticle trajectories with different numbers of An- 
dreev reflections. The magnitude of SGS decreases both 
at T ^ and T ^ Tc, which resembles the behavior of 
SGS in the I-V characteristic of long ballistic SNS V^c- 
tion with perfect interfaces within the OTBK modelll3 A 
small "residual" SGS in cuxrent noise, similar to the one 
in the TV characteristicjlj should occur at T ^ due 
to normal scattering at the interface or due to proximity 
effect [see comments to Eq. (p^)]. 
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FIG. 2. Spectral density S of current noise vs voltage (a) 
and its derivative dS/dV vs inverse voltage (b) at different 
temperaifcares. Dashed line shows the result for normal metal 
junctionkl at T = 0.3Tc. 

The role of inelastic scattering is to suppress the 
nonequilibrium of the subgap electrons created, by 
MAR. Within the relaxation time approximation,Blla the 
nonequilibrium distribution of the subgap quasiparticles 
holds as soon as the inelastic relaxation time Te{E) at 
the characteristic energies i? ^ A is larger than the diffu- 
sion time through the junction, Td{V) ~ {2A/eVYd^ / D. 
Thus, in this model, the enhancecL-S|hot noise (as well as 
the SGS in the I-V characteristicEj) can be observed if 

— 1/2 

the apphed voltage is sufficiently large, eV > 2AWe , 
where = E'y\,t^[A) /Ti\ at smaller voltages, S{V) 
should decrease and approach the thermal noise level. 
The noise temperature is equal to the physical tempera- 
ture T if the inelastic scattering is dominated by electron- 
phonon interaction (assuming that the phonons are in 
equilibrium with the electron reservoir). If the electron- 
electron (e-e) scattering dominates, the noise tempera- 
ture may exceed the temperature T of the electron reser- 
voir if this temperature is small, T -C A (hot electron 
regime). The reason is that at low temperature, the sub- 
gap electrons are well decoupled from the reservoir (elec- 
trons outside the gap) due to weak energy flow through 
the gap edges. 

In order to quantitatively analyze this situation, we 
consider the small voltage limit, eV <^ A, in equation 
Eq. (^), taking into account the e-e collision term. Due 
to weak spatial dependence of the population numbers 
at small voltage, they can be replaced by their boundary 
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values, rig' {E) « (E) = n{E), in the e-e collision 
integral lee (n) ■ This allows us to easily include the colli- 
sion term in the recurrences of Eq. . Within the same 
approximation, these recurrences are to be considered as 
differential relations, which results in the diffusion equa- 
tion for n{E), 



dE^ 



Iee{n). 



(15) 



where De — (eT^)^i?Th/?i is the effective diffusion coeffi- 
cient in energy space. The finite resistance Rns of the NS 
interfaces, which partially blocks quasiparticle diffusion, 
can be taken into consideration by renormalization of the 
diffusion coefficient, De De[1 + (d/^o)(i?ws/i?)^]"^ 
(see Ref. |l|). 

Equation (|l^) describes the crossover from the "colli- 
sionless" MAR regime to the hot electron regime as func- 
tion of the parameter Z)£;ree(A)/A^. In the hot electron 
limit, 3> DETee, the collision integral dominates in 



Eq. (15), and therefore the approximate solution of the 
diffusion equation is the Fermi function with a certain 
effective temperature Tq ^ A. The value of Tq can be 
found from Eq. ( |l5|) integrated over energy within the 
interval (—A, A) with the weight E, taking into account 
the boundary conditions n(±A) = nF(±A) and neglect- 
ing the exponentially small derivative dn/dE at the gap 
edges. At zero temperature of the reservoir, we obtain 
an asymptotic equation for Tq, 



{eVfW, exp(A/To) = roA(l + Tq/A), 



(16) 



which shows that the effective temperature of the subgap 
electrons decreases logarithmically with decreasing volt- 
age. The noise of the hot subgap electrons is given by 
the Nyquist formula with temperature Tq , 



S{V) = (4ro/i?) [1 - 2exp(-A/To) 



(17) 



where the last term is due to the finite energy interval 
available for the hot electrons, \E\ < A. Equations (16), 
( p^ ) give a reasonably good approxiraation to the result 
of the numerical solution of Eq. ( [l5| ) .til 

In summary, we have calculated current noise in a 
long diffusive SNS structure with low-resistive interfaces 
at arbitrary temperatures. Whereas the I-V charac- 
teristic is approximately described by Ohm's law, the 
current noise reveals all characteristic features of the 
MAR regime: "giant" enhancement at low voltages, pro- 
nounced SGS, and excess noise at large voltages. In the 
limit of strong electron-electron scattering, the junction 
undergoes crossover to the hot electron regime, with the 
effective temperature of the subgap electrons decreasing 
logarithmically with the voltage. 
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